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The unifying theory of scaling in thermal convection (Grossmann & Lohse (2000)) sug- 
gests that there are no pure power laws for the Nusselt and Reynolds numbers as func- 
tion of the Rayleigh and Prandtl numbers in the experimentally accessible parameter 



regime. In |Grossmann fc Lohse (2001) the four dimensionless parameters of the the- 
ory were fitted to 155 experimental data points by Ahlers & Xu (2001) in the regime 
3 x 10 T < Ra < 3 x 10 9 and 4 < Pr < 34. Meanwhile the theory is on one hand well con- 
firmed through various new experiments and numerical simulations. On the other hand 
these new data points provide the basis for an updated fit in a much larger parameter 
space. Here we pick four well established (and sufficiently distant) data points Nu(Ra,Pr) 
to obtain the four dimensionless constants q of Grossmann and Lohse's unifying theory 
and show that the resulting Nu(Ra,Pr) function is in agreement with almost all estab- 
lished experimental and numerical data up to the ultimate regime of thermal convection, 
whose onset also follows from the theory. In addition, Re(Ra,Pr) is provided. 



1. Introduction 

Thermal convection is omnipresent in science and technology and its paradigmatical 
representation is Rayleigh-Benard (RB) convection: a fluid in a sample heated from be- 
low and cooled from above. This system has received considerable attention in the last 
decades ( | Ahlers e£~aLl|20096| [Siggia 1994] |Lohse fc Xia||2010| ), with one focus on the 
scaling properties of the global heat transport of the system. A now widely accepted 



viewpoint is the Grossmann-Lohse (GL) theory (Grossmann fc LohsepOOO, 2001, 2002 



2004). The basis for this theory of scaling in RB convection are exact global balances 



for the energy and thermal dissipation rates derived from the Boussinesq equations and 
the decomposition of the flow in boundary layer (BL) and bulk contributions. The scal- 
ing of the dissipation rates in the BLs is assumed to obey Prandtl-Blasius-Pohlhausen 
scaling (Schlichting 1979), which is justified as long as the shear-Reynolds numbers of 
the BLs are not too large, and the scaling relations in the bulk are estimated based 
on Kolmogorov-type arguments for homogeneous isotropic turbulence. While the theory 
gives the different scaling relations for the individual contributions to the energy dissi- 
pation rates in the bulk and in the BL, namely e u ^ u lk an d £ u ,bl, and to the thermal 



dissipation rates in the bulk (background) and in the BLs (plus the plumes, see Gross- 
mann & Lohse (2004)), namely eg^ u ik and cq^l, the absolute sizes of these four relative 
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contributions are not given by the theory. They are expressed in four dimensionless pref- 
actors Ci, z = 1,2,3,4 for c UjB l, e Ul buik, £o,bl, and e e ,buik, respectively, which have to be 
adopted to experimental or numerical data for Nu(Ra, Pr). 

When the theory was developed early this century, such data were scarce and often 
contradicting each other, due to sidewall and plate effects, insufficient knowledge of the 



material properties of the fluid, lack of numerical resolution and other problems. Gross- 



mann & Lohse 

3 x 10 Y <la<Tx 10 9 and 4 



p00l| used 155 data points for Nu(Ra, Pr) in the pa r amete r range 
< Pr < 34 obtained by |Ahlers fc Xu| ( |2QQl| ), which 
was the most extensive data set at that time. This fixed Nu(Ra, Pr) for all Ra and Pr, 
considered as valid up to the meanwhile found (He et al. ( |2012o| )) ultimate regime of 
thermal convection, where the Prandtl-Blasius type BL becomes unstable. Re(Ra, Pr) 
was fixed (cf. Grossmann fc Lohse] ( |2QQ2| ) with one extra adoption of the prefact or a in 
the Prandtl-Blasius scaling relation A M = aL/y/Re to the experimental data of |Qiu fc] 



Tong (20011), where A n is the mean thickness of the kinetic BL and L the height of the 



sample. 

Although the data to which we adopted the four prefactors q were relatively local 
in parameter space, the theory was rather successful in describing the global behavior 
Nu(Ra, Pr) and also Re(Ra, Pr), as described in detail in Ahlers et al. (20096). This 
included the prediction that for Pr « 1 the onset to the ultimate regime should take 
place when Ra is of the order of 10 14 . This prediction was based on an assumed onset 
of a sheared BL instability at a shear Reynolds number Re s « 420, which is the value 
given in |Landau fc Lifshitz (1987). Indeed, very recently He et al (20126) have found 
the onset of the ultimate regime at this very Rayleigh number. 

Thanks to joint efforts of the community the experimental and numerical data situa- 
tion for Nu(Ra, Pr) has considerably improved in the last decade. Measurements have 
been extended to a much larger domain in the Ra-Pr parameter space, see the updated 
phase diagrams in figure [I] and figure [7] and plate- and sidewall corrections are much bet- 
ter understood and taken into account (|Brown et al] (|2005|); lAhlersI (|2000|); iRoche et al 



(2001);|Verzicco|(|2002|);|Niemela & Sreenivasan (2003); Ahlers et al 



( |2009fr| )). Further- 



more, due to the increasing computational power and better codes the numerical data are 
now well converged, confirming and complementing the experimental data. Meanwhile 

achieved Ra = 2 • 10 12 at Pr = 0.7inar = l/2 sample and 



Stevens et al (2010c 



2011a 



obtained a good agreement with the experimental data of |He et al\ ( |20126| ) and |Niemela 



et al (2000). 



This situation calls for a refit of the four prefactors q of the GL theory, in spite of the 
success of the theory with the prefactors of Grossmann & Lohse ( |2001 ): It is clear that 
the surface Nu(Ra, Pr) above the Ra-Pr parameter space will be much more stable and 
"wobble" less if we put it on four distant and trust able "legs" Niti(Rai, Pr^), i = 1, 2, 3, 4, 
rather than putting it on four "legs" somewhere in the center but close to each other. 
We emphasize that we obviously only need four such "legs" for the four constants. The 
shortcoming of the old set of q was particularly obvious for small Pr, say Pr < 1 (see 
figure [5]), because at the days of |Grossmann fc Lohse (2001) no reliable information was 
available in that parameter regime and therefore no Nusselt data of that regime had been 
included into the fit. 

The structure of the paper is as follows: In section 2 we will provide the refit of the GL 
theory for an aspect ratio T = 1, leading to Nu(Ra, Pr) in the whole parameter space 
up to the ultimate state. In section 3 we discuss the robustness of the fit. In section 3 we 
will show that this fit also describes the available data for T = 1/2 and will in particular 
discuss the onset of the ultimate regime. Section 4 gives conclusions and an outlook of 
the new challenges. 
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Ahlers and Xu (2001) 

Xia et al. (2002) 

Fleischer and Goldstein (2002) 

Niemela and Sreenivasan (2003) 

Sun et al. (2005) 

Funfschilling et al. (2005) 

Roche et al. (2010,r=1 .14) 

Burnishev et al. (2010) 

He et al. (2012) 

Emran and Schumacher (2012) 
Stevens et al. (2010-2013) 



2001, 2002 



lane for RB convection according to Grossmann and 
2004)) in a Y — 1 sample with no-slip boundary 
; the lower nearly parallel solid line corresponds 
to, bulk] and along the long-dashed 



Figur e 1. Phase diagram in Ra — Pr 
Lohse (Grossmann & Lohse (2000 
conditions. The upper solid line means He 
to € u ,bl = £u,buik\ the curved solid and dashed line is 6q,bl 

line X u = A,?, i.e., 2aNu = y/Re. The dash-dotted line indicates where the laminar kinetic BL 
is expected to become turbulent, based o n a critical s hear Rey nolds number Re* — 284 of the 

~al\ ([T997t; [Glacier et al\ ([19991; | Ahlers & 



Xu 



kinetic BL, see text. The data are from |Cioni et 

(12001|); |Fleischer fe Goldstein| (|2002|) ; fxia et al.\ (120020; l^haumat et al.\ ^2QQ2tt |fJiemelaT 



Sreenivasan! ( 12003D;|Sun et a/.|Q2005l );|Funfschilling et al (2005) ;TRoche et a/.|l20T0|); |Burnisnev 
et a/.|Q2010t;|Kmran & Schumacher] q2012p ; |He et al .| ( |2012at ; |Stevens et ai\ ( |20106|al ); |van"c[er 
Poel et al\ ( |2013| ). Note that lor the Stevens et al. data points from different papers nave been 
combined in the graph. 




Figure 2. Compensated three-dimensional visualization of a) Nu(Ra,Pr) and b) Re(Ra,Pr). 
The four points used to fit the GL parameters c\ — 112.3161, C2 = 67.6078, C3 = 0.9318 and 
C4 = 0.0921 have been indicated by the black points in the Nu(Ra,Pr) graph. 
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Figure 3. Comparison of the i?a-scaling of the original GL-fit from |Grossmann &; Lohse (2001 ) 
(a) with the new fit (b) for water, i.e. P r = 4.38 unles s mentioned other wise, in a 1' = 1 sample 



The circles dFunf schilling et aZ.N'2005)) and squares ( |Sun &; Xia| (|2005)) indicate experimental 
results. Open symbols indicate the uncorrecte d data and solid symbol s the data after correction 
for the fin ite plate conductivity. The diamonds jLakkaraju et a l. (2012), Pr — 5.4), rig ht pointing 
triangles (Stevens et al. (20116)), and left pointing triangles ( |Shishkina & The ss (2009)) indicate 
results from numerical simulations. 

2. Refit of the GL theory for r = 1 

The GL theory describes Nu(Ra,Pr) and Re(Ra,Pr) with the following two coupled 
equations flAhlers et~al\ ( |2009fr| ) ) , 

Re 2 



(Nu - l)RaPr~ 2 = a 



g(^/Re L /Re) 
Nu — 1 = c 3 Re 1/2 Pr 1/2 { f 

c 4 PrRef 




(2.1) 



(2.2) 



where the crossover functions f and g model the crossover from the thermal boundary 
layer nested in the kinetic one towards the inverse situation and from A n = aLj\fRe 
toward A M ~ L, respectively; for details, see Grossmann & Lohse (2001). As described by 



Grossmann & Lohse (2002) the prefactor a = 0.482 is obtained from the experimental 
data of Qiu & Tong (2001) and in addition we use ReL = 1.0, where ReL is that Reynolds 



number for which the BL thickness according to Prandtl-Blasius is of order of the sample 
extension. 

In order to get accurate values for the four dimensionless prefactors q it is necessary to 
provide four data points with as much information about the richness of the RB system 
as possible, which means that data points from different regimes should be selected. 
Therefore we determined the q from the data points of Funfschilling et al. ( 2005) at 
Ra = 1.8 x 10 7 and Ra = 2.25 x 10 10 , both with Pr = 4.38, the dat a point fromjXia et al. 
( 2002[ ) with Pr = 818 at Ra = 2.04 x 10 8 , and the data point from |Cioni et al] fll997| ) at 
Ra = 1 x 10 T with Pr = 0.025. The location of these data points in the RB phase diagram 
is indicated by the large red squares in figure [l] and by the black dots in the corresponding 
three-dimensional Nu(Ra,Pr) visualization in figure [2^t. Figure [l] shows that these data 
are indeed within different regimes. The reason for choosing these specific data points is 
two- fold. First of all we consider these four data point to be reliable. And apart from the 
data point by Xia et al. (2002), which is the only experiment in that large Pr regime, all 
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Figure 4. Comparison of the i?a-scaling of the original GL-fit from Grossmann & Lohse (2001) 
(a) with the new fit (b) in a F = 1/2 sample and varying Pr, see ph ase diagra m in figure 
u\ The right p ointing triang l es are the experimental data from 



Castaing et q/.| d!989> with 

wall corrections |Roche et al.\ (|2010|), left pointi ng triangles ( Fleischer & Goldstein| (|2002|)), up- 



ward p ointing triangles (Chaumat et al. (2002)), do wnward p ointing triangles (Chavanne et al 
'2001)), diamonds (Nie mela et al.\ (|2000|)), crosses (jAhlers et al. (2009a)) , stars ( |Hoche et cu 
2010)), circles |He et a/.|(|20126D;|Arilers et j20126|) ), and plusses ) Urban et a/.|(|20111|2012 



indicate experimental data and the squares results from numerical simulations (Stevens et al. 
( 2010cl[20lT^ )). 




Figure 5. Comparison of the Pr-scaling of the original GL-fit from I Grossmann fc Lohse|( 2001) 
(a) with the new f it (b) for different Ra . The right pointing triangle ( |Cioni et al.\ (|1997[)) , left 
pointing triangle jHoranyi et al\ (|1999|), and the upward pointing triangle ( jHossby 



(1969 



indicate e xperim ental results. The circles ( |van der Poel et al. \ ( |2013| ) and squares ( |Verzicco &| 
Camussi| (|1999|)) indicate numerical results obtained in a cylinder with aspect ratio F — 1. The 



[Ca^ 

plus: 



'I 



plusses ( jPetschel e t al. (2013)) indicate numeric al results obtained i n a periodic domain and the 
diamonds ( |Kerr & Herrin g (2000)) and crosses ( |Breuer et al. (2004)) numerical results obtained 
in a box with free slip boundary condition at the sidewall. The colors black, purple, magenta, 
red, and blue (from bottom to top) corresponds to the Ra numbers 5 x 10 5 , 10 6 , 5 x 10 6 , 10 7 , 
and 10 8 . 



points agree very well with experimental or numerical data from other groups, see figures 
[3]and|5] In addition, these four data points are relatively far apart in the Ra-Pr parameter 
space to ensure that they provide the theory with as much information about the richness 
of the RB physics as possible. To provide infor mation about the i?a-scaling we selected 
the measurements of Funfschilling et al. ( 2005| ) at Ra = 1.8 x 10 7 and Ra = 2.25 x 10 10 
with Pr = 4.38. In order to include information about the transition between the 'upper' 
and 'lower' regimes, which is modeled by the crossover functions / and g, it is necessary 
to provide data points in the low, intermediate, and high Pr number regime. We do this 
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selecting next to the intermediate Pr n umber data from Funfschilling et al. (2005), the 
low Pr = 0.025 number meas urement by|Cioni et al. (1997) at Ra = 1 x 10 7 and the high 
Pr = 1352 measurement by |Xia et al\ (|2002D at~Ra = 1.78 x 10 9 . Altogether the four 



data points provide information from three different Pr numbers and four different Ra 
numbers. From these four data points we determine the q with a fourth order Newton- 
Raphson root finding method and also by using a trust-region-reflective optimization. 

Both methods give c x = 112.3161, c 2 = 67.6078, c 3 = 0. 9318 and c 4 = 0.0921 

In figures [3] to [5] we compare the GL-fit determined by Grossmann & Lohse (2001 ) with 
this new GL-fit, These figures clearly reveal that the new GL-fit is much closer to the 
data in the low Pr number regime, while maintaining the similar excellent agreement for 
the high Pr number data as before. We emphasize that this excellent agreement with the 
data from experiments and simulations that are not included in the fitting procedure, 
i.e. only four data points are used in the fitting procedure, confirms that the values 
we found describes Nu(Ra, Pr) well in the regime that is nowadays covered by state of 
the art experiments and simulations. It is also noteworthy that figure [3] and [4] show that 
the Ra number scaling is perfectly predicted by the GL-theory for Ra values that are 
decades higher than the highest Ra number point that is used to determine the q values, 
i.e. Ra = 2.25 x 10 10 , thus showing the predictive power of the GL-theory. 



3. Robustness 

To illustrate the robustness of the fit presented above, we show the results of a fit 
through four other data points, namely the data points from | Funfschilling et al. ([2 005) 



at Ra = 2.96 x 10 7 and Ra = 1.92 x 10 1U with Pr = 4.38, the one from Xia et al 



2002 



at Ra = 2.24 x 10 with Pr = 554 and finally the data point by|Kerr fc Herring 



2000[ ) at Ra = 10 7 with Pr = 0.07. Three out of these four data points lie relatively 
close to the original four data points, but the low Pr = 0.07 point from |Kerr fc Her ring 
(2000) substantially differs from the original Pr = 0.025. The reason that three of the 
four points are close to the original four points in the Ra — Pr parameter space is that 
one can only select "reliable legs" in regimes were many measurements have been done 
and these regimes only cover a limited part of the parameter space. 

The resulting GL coefficients are c x = 114.1135, c 2 = 38.0299, c 3 = 0.9226 and c 4 = 
0.0677, compared to c t = 112.3161, c 2 = 67.6078, c 3 = 0.9318 and c 4 = 0.0921 of the 
fit described above. In order to compare the two fits we compare the relative difference 
in Nu(Ra, Pr) calculated in the fit described in the previous section and Nu calculated 
from this additional fit in the parts of the parameter space where the GL fit is valid. A 
comparison between both fits shows that the difference is very minor in the regimes IV U , 
II U1 and 7 n , and that the differences increase in the regimes 11^ IVi, and III U , which 
are very far away from the region in the parameter space where reliable data points are 
available. The reason is that a very small variation in the measurements point can lead 
to significant differences if the implied information is extrapolated over many decades in 
Ra and Pr using the GL-theory. For the fits compared here the differences increase up 
to about 10%. 



4. GL theory for T = 1/2 and ultimate regime 

In principle, the q depend on the aspect ratio T. However, it is well known that only 
small differences in Nu are observed between T = 1/2 and T = 1 (Ahlers et al. (20096)). 
This weak aspect ratio dependence is confirmed by figure [4j which shows that the Ra 
number scaling for Pr = 0.7 in a T = 1/2 sample is captured very accurately by the new 
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Figure 6. Relative difference between Nu calculated from the original fit and Nu calculated 
from the additional fit. The color scale ranges from blue to red, signifying 0% to 10% difference 
respectively. 



fit for T = 1, and in the low Ra number regime the new fit is even much better than by 
the original fit of Grossmann & Lohse (2001). 

The location in Ra-Pr space of the various regimes of the GL theory is based on the 
coefficients q. The updated lines that encompass the regimes are plotted in the phase 
diagrams shown in figures [T] and [7] The line that indicates the onset of the ultimate 
regime, where the kinetic boundary layer has become turbulent, i s now b ased on the new 
Ci and the transition at Ra = 5 • 10 14 , observed by He et al. (20126) for Pr = 0.86. 
This gives i?e* = 284 instead of the previously used i?e* = 420 taken from pipe flow 
(Landau & Lifshitz (1987)). This new value for Re* s is compatible with observations for 
flows along plates by Hansenj ( |1928 ) and with the value 320 promoted by Niemela & 



Sreenivasan (2003). 



The phase diagram in figure [7] shows that the measurements of He et al. ( 2012fr| ) up 
to Ra ~ 10 15 at Pr = 0.86 are the only experiments that have reached the ultimate 
regime. They observe the onset of the ultimate regime at Ra = 5 • 10 14 and a transition 
region for 10 13 < Ra < 5 • 10 14 . The experiments by He et al. (20126) are the only room 
temperature experiments for Ra > 10 12 , while all other experiments that have reached 
these Ra numbers are low temperature experiments with Helium close to the critical point 



(Chavanne et al. (1997, 2001 


); 


Niemela et al.\( 


2000 2001 


); 


Niemela & Sreenivasan ( 2006 ) ; 


Roche et a/.|fl2010|); Urban et al. (2011 2012 


)). In these low temperature experiments it 



is difficult to reach the ultimate regime because the Pr number increases with increasing 
Ra, see figure [7] Nevertheless the low temperature experiments by Niemela et al. (2000) 



seem to come very close to the ultimate regime and one may wonder why the transition 



region observed by |He et al. (20126) was not observed in the Niemela et al. (2000) 



experiments. We believe that the scatter of the |Niemela et al ( 2000| ) data at this highest 
Ra due to the uncertainties in the fluid properties in combination with the fact that the 
transition is smooth is the reason for this. Figure [4] shows that the scatter in the Niemela 



et al. (2000) data is similar to the trend observed in the transition regime. The phase 



diagram shows also shows that other low temperature experiments by |Chavanne et al. 



( 1997||200lD , |Roche et a/.| ( |2010D , and |Urban et a/.| ( |2011[|2012D do not reach the ultimate 
regime and therefore no transition to the ultimate regime due to a BL shear instability 
is expected in these experiments. 
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Castaing et al. (1989) 
Chavanne et al. (1997) 
Glazier et al. (1999) 
Niemela et al. (2000) 
Ahlers and Xu (2001) 
Chaumat et al. (2002) 
Burnishev et al. (2010) 
Ahlers et al. (2009) 
Roche et al. (2010) 
Stevens et al. (2010-2011) 
Urban et al. (2011-2012) 
He et al. (2012) 



Figure 7. Phase diagram in Ra — Pr plane for RB convection in a Y — 1/2 sample with no-slip 



(1989); Chavanne et a/.|fll997|);|Glazier et al. 


(1999); Niemela et al (2000 


);|Ahlers & Xufi 


20011); 


Chaumat et at. 


(2002); Ahlers et al (20096); 


Burnishev et al 


(2010); Roc 


tie et al\ (12010); 


Urban 


et ai. (|2011||2012|);|He et ai. 


(20126|);|Steven 


s et al 


(|2010c||^ 


>0lTaf. 



5. Conclusions and outlook 

In this paper we have used the availability of new experimental and numerical data, and 
our increased understanding of the physics of the Rayleigh-Benard system to determine 
the prefactors of the unifying theory for scaling in thermal convection, i.e. the Grossmann- 
Lohse theory, much more accurately. The resulting Nu(Ra,Pr) function is in very good 
agreement with almost all established experimental and numerical data up to the ultimate 
regime of thermal convection, and has significantly improved the predictions. In figure [4] 
one can notice the onset of the ultimate regime in the Nu(Ra) scaling of the measurements 
of He et al. (|20126| ). Extensions of the GL theory to the ultimate regime (see |Grossmann| 
& Lohse ( 2011[ )) are able to explain the observed Reynolds number scaling in the ultimate 
regime as well as the origin of the log-profiles Grossmann & Lohse ( 2012| ), observed in 



the ultimate regime Ahlers et al. (2012a) 

In line with Grossmann & Lohse (2001), we have determined the prefactors from 



experimental measurements. This has great value as it shows that the information of 
only four data-points is sufficient to accurately predict Nu(Ra,Pr) up to the ultimate 
regime. All is based on the GL theory, which builds on exact global balances for the 
energy and thermal dissipation rates, derived from the Boussinesq equations, and the 
decomposition of the flow in boundary layer and bulk contributions. 

A further challenge we want to pursue is to calculate the q directly from the fluid equa- 
tions, without the input of any experimental or numerical data, or at least quantitatively 
relate their values to important fluid concepts like Prandtl-Blasius-Pohlhausen theory, 
the von Karman-Prandtl theory, etc. in order to get an even deeper understanding of the 
GL theory. 
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